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ABSTRACT 


For heavy nuclei, very few measurements of the nuclear giant multipole 
resonances have been conducted using inelastic electron scattering. In 
this work, measurements were made to identify the giant multipole resonances 
of ae Five inelastic scattering spectra were obtained and analyzed for 
beam energies and angles of 75 MeV / 75. , 75 Mev / 90°, 75 Mev / 120m 
90 MeV / 75', and 90 MeV / 90°. In each spectrum, giant multipole rasotces 
were observed at excitation energies of 7.1, 9.1, 10.9, and 14 MeV. The 
strengths and multipolarities were found to be B(M2) = 15,000e2fm*, 

B(E2) = 8,000e2fm4, B(E2) = 9,000ce2fm4, and B(El) = 90e2fm2, respectively. 
The full widths at half maximum for these resonances were found to be, 
pespecuively. G2le07-9, 3.1, andg459 MeV. 

The resonances found at 10.9 and 14 MeV confirm the only other measure- 
ments made on gold | 4]. The resonance found at 7.1 MeV agrees with those 
findings concerning the magnetic character of the transition but differs in 
the multipolarity Aas. WEB. our assignment may also be consid- 
ered tentative since insufficient data were obtained at the backward scat- 
tering angles. The resonance found at 9.1 MeV has not been previously 


reported. 
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ieee! RODUCTION 


The giant dipole (El) resonance has been well known for some time. 
Recent experiments using inelastic electron scattering have revealed the 
existence of other giant resonances in a number of nuclei and in this paper 
we report our investigations about these resonances. 

Many inelastic scattering measurements have been made in the past 15 
years; however, the giantresonances considered in this paper appear as small 
peaks superimposed on a large continuous spectrum of electrons called the 
elastic radiation tail. Low background and high accelerator beam stabil- 
ity are required to detect these resonances and it is only recently that 
they have been investigated. Studies using Ce, La, Pr [ Refs, 1,2] and 
Pb [ Ref. 3] indicate that heavy nuclei tend to have an El resonance at 
about 13 to 15 MeV and an E2 resonance at about 10 to 12 MeV. Other reso- 
nances, some of which were magnetic in character, were observed at excita- 
tion energies ranging from 7 to 9 MeV. Unpublished work at Darmstadt 
[ Ref. 4] on 177 au has also confirmed this trend: an El resonance was found 
at 10.8 MeV; and one other resonance was located at 7.1 MeV and tentatively 
identified as Ml. Another: recent experiment on 907, [ Ref. 5] found a giant 
resonance (EO or E2) at 28 MeV excitation in addition to the usual E2 and 
El transitions near 15 MeV. 

The investigations reported here were the first attempt to measure 
glant resonances at the NPS LINAC. The procedures and data handling tech- 
niques necessary for the study of giant multipole resonances by inelastic 
electron scatterings have been developed and refined. 

Experimental values for the inelastic scattering form factors were 


obtained using incident energies of 75 and 90 MeV and scattering angles of 





75, 90, 120, and 150 degrees. The results of this work agree with other 
reported works [ Ref. 4] in that the previously reported resonances were 
observed. A previously unreported resonance was found at 9.1 MeV with 
B(E2) = 8,000e2fm*. Additionally, preliminary indications for the exis- 
tence of an E2 resonance in the vicinity of 23 MeV were obtained. 

This thesis is organized into five sections and one appendix. In 
Section II a brief review of inelastic scattering and giant resonance 
theory is presented. Section III briefly describes the NPS LINAC and the 
experimental procedure used by the authors. Section IV discusses the meth- 
ods employed for data reduction. Topics include radiation corrections, 
calculations of the elastic radiation tails, procedures used in the loca- 
tion and characterization of the giant resonances, and error analysis. 
Experimental results and conclusions are presented in Section V. The cal- 
culations involved in this work are quite complex and cannot reasonably be 
performed without the use of various sohpisticated computer programs. The 
appendix contains the program listings and detailed explanations concerning 


their use. 





tie LHEORY 


A. INELASTIC ELECTRON SCATTERING 

When a relativistic, mono-energetic electron beam is incident on a 
thin target most electrons pass through, but some collide elastically with 
target nuclei and a few collide inelastically, thereby exciting nuclei to 
higher allowed energy states. 

The energy distribution of the scattered electrons is determined ex- 
perimentally at various scattering angles. Most of the electrons observed 


at a given energy, E less than the incident energy, Eas are due to elas- 


£? 
tic collisions with nuclei accompanied by energy loss in passing through 
the target. Energy loss mechanisms include: (1) ionization losses, 
(2) emission of a photon during the scattering event, and (3) bremsstrah- 
lung while passing through the target both before and after the elastic 
scattering event. This "elastic tail," and any other "background" affect- 
ing the observed spectrum, must be calculated and subtracted from the ex- 
perimental energy spectrum to obtain the inelastic spectrum. The resulting 
unfolded inelastic spectrum consists of the inelastic peaks and their 
associated tails. Each inelastic peak corresponds to an excited state of 
the uel ue: 

Inelastic collisions are classified into three types, according to the 
mechanism of interaction with the nucleus: 

1. Surface Nucleon Excitation 

One or more surface nucleons are excited with excitation energies 

ranging from several KeV to a few MeV. The higher energies tend to be 
associated with collective excitations and have larger characteristic 


resonance widths. 
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2. Giant Resonances 
Many nucleons participate in oscillations causing the giant reso- 
nances. These resonances have very large natural widths and are found 
typically in the region of ten to 20 Mev excitation forvheavy nucHei. This 
large resonance width occurs because the excitation energy is large enough 
that the excited nucleus can decay by particle emission. The short life- 
times associated with particle emission result in the large widths because 
of the Heisenberg uncertainty relationship between energy and time. 
3. Quasi-Elastic Resonances 
In this process a nucleon is ejected from the nucleus during the 
collision. The incident electron interacts with essentially only one 


nucleon. These 'quasi-elastic'' collisions produce a very wide peak in the 


spectrum at an excitation energy of a few hundred MeV for heavy nuclei. 


B. GIANT RESONANCES 

According to the Generalized Goldhaber-Teller Model | Ref. 6]. nuclei 
are capable of "collective oscillations" in which all the nucleons oscil- 
late inside the nucleus. There are four separate fluids in motion: neutrons 
with spin up, protons with spin up, neutrons with spin down, and protons 
with spin down. Each fluid is characterized by the spin and the isospin 
of nee particles. 

The motion of pie (egue fluids, two of which are 180° out of phase with 
the other two, can be characterized by two numbers, S and T. If protons 
and neutrons move in phase, T=0; if they move 180° out of phase, T=1l. If 
particles with spin up are in phase with those with spin down, 5=0; if they 
are 180° out of phase, S=l. A transition to an S=1 state changes the mag- 
netization and is therefore called a magnetic transition. A transition to 


an S=0 state is called an electric transition. 
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The electric and magnetic fields of nuclear excited states are described 
in terms of multipole expansions. A transition is labeled according to the 
term in the multipole expansion which is predominantly affected. For ex- 
ample, a transition from the ground state meetin: donee the magnetic (or 
Piletric) dipole term in the electromagnetic multipole expansion is labeled 
M1 (or El). If the electric (or magnetic) quadrupole term is involved, the 
transition is labeled E2 (or M2). An electric monopole transition is labeled 
EO. There are no magnetic monopole transitions. In general, L units of 


angular momentum are required for an EL or ML transition from the ground 


state. 


C. NUCLEAR FORM FACTORS 
1. Development of the Form Factor 
In the plane wave solution to the non-relativistic wave equation 
we find that the differential elastic scattering cross section can be ex- 
pressed as the magnitude of the square of some function of scattering angle, 
£(@), where 
f(0)=- qf UCP) ead i. 

a oe 
and U(r) = 42 V(r), where V(r) is the potential of the electron in the 
field of the nucleus, and K is the wave vector of the incident wave. 
Replacing U(r) by the related electrostatic potential, d(r), through the 


relation V(r) = q h(x), we have 
f(a) =+2me fre’ plFrd?r. 


We can relate @(r) to the charge distribution p(r) through Poisson's equa- 
tion; taking the fourier transforms and equating integrands leads to the 


result 


EZ 





f(o)= BME (2m) 4% vf), 


with 
ple): (a) ree dr = (ae) ZeF(e), 


Hence, 


‘_ 2me’Z 
f (9) = RR le) . 


We can now solve for the’differential cross section to obtain 


do _ amee re 
oe. = FE. 


For a point charge it can be shown that giz Amey , so we have 


do = Clo K 4 
de I, Fel 


We call F(K) the form factor for the charge distribution with 


F(k) = ze [ome ars a (II-1) 
For inelastic scattering, the charge distribution changes during the 
transition. Inthis case, “p(r) is replaced by the four-vector Dirac Are 
rent operator, Jus associated with the transition involved. 
Jig Mott Cross Section 
The differential cross section per unit solid angle for elastic 
scattering of a relativistic electron by a point nucleus of charge 2 was 


shown by Mott [Ref. 7] to be: 


sal z?e*cos"6/2 /R= = (60 
do}, 4E; sin*e/2 dQ] mott 


where R is the reccil factor defined by the relation 
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L 


1+ 2E; 
*Mc?sin ? 6/2 


R 


Using the concept of the form factor we are able to generalize this result 


to finite charge distributions to obtain 


Gal (dal oe 


For heavy nuclei where E;< 100 MeV, R=13; therefore 


Fal - aA mot a | 


If this expression is used with the form factor given by equation II-l, 


F®)|*R. 





(II-2) 





we have the Born approximation for elastic scattering. For large Z and low 
energy this approximation is inaccurate and the cross section must be cal- 
culated by partial wave analysis [ Ref. 8]. 
3. Electron Scattering from a Proton 
The electromagnetic interaction Hamiltonian between the electron 


and a point proton is expressed in four-vector notation as He .= ~ dpa. 


Calculating Ju in terms of the Dirac theory operators [Ref. 9| pives 
Hy.= 2 ie (Mixx) (ie 72) (HK Te) 
ji 


Solving the Dirac wave equation yields the point Dirac proton solution: 


do _ [de _g_ do 
do - Gal mort R (+ ginger tong?) = {ain 


The first term in the parentheses corresponds to the electric interaction 





| mott 


and the remaining term corresponds to the magnetic interaction. General- 
a 
izing to finite charge distributions and adding another term, Fy(q°), 


account for the anomalous magnetic moment distribution we obtain 
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9a =fa5) moni (Isl + a am?c? 


MncmElrst two Cermseaccount for the electric interactions#@and the last term 


aR 4 a amic7tan’4/2 | F +B), 


4am?c 





accounts for the magnetic scattering; q? is the square of the four-momentum 
transfer and “is the proton magnetic moment in nuclear magnetecns. 

This result, although it does not apply directly to these investigations, 
does show the presence of two form factors, due to scattering from the charge 
and magnetic moment of the proton. 

4. Plane Wave Born Approximation (PWBA) 

The Born approximation is used to calculate inelastic scattering 
cross sections based on an assumed nuclear transition charge distribution. 
The PWBA uses first order perturbation theory to calculate cross sections 
for interactions between an electron and the electromagnetic field of the 
nucleus. A multipole expansion is used to describe the fields. The tar- 
get nuclei are assumed to recoil non-relativistically. The PWBA assumes 
plane wave solutions for the incident and scattered electron. 

For E;, Eg>>Ex, Isabelle [Ref. 10] has shown that the differential 
inelastic cross section per unit solid angle is given by 


fae * ide) y Ret +(3 +tan?6/2) Fiq,L)| | 


and 


do do a 2 2 
(Sal = (ge Revst ( > tfan 6/2) Fab) 


where Fp and F, are the form factors for longitudinal and transverse photon 


L 
absorption. For elastic scattering, the PWBA results are equivalent to 


equations II-l and II-2. 
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5. Distorted Wave Born Approximation (DWBA) 


For heavy nuclei, the coulomb field surrounding the nucleus dis- 
torts the electron wave function sufficiently that a plane wave solution 
is inadequate. The DWBA assumes phase~shifted spherical waves for the in-~ 
cident and scattered electron wave functions. Unlike the plane wave Born 
approximation, a form for the charge distribution function must be assumed 
before calculations are begun. 

For elastic scattering, the cross section may be calculated using a 
Fermi charge distribution. However, this is not appropriate for inelastic 
scattering, in which the outer nucleons only would be expected to take part 
in a nuclear transition. A good approximation was suggested by Tassie 
[Ref. iy in which the transition charge distribution is assumed to be 
the derivative of the Fermi function. Computer programs for calculation 


Sreeiewaitierential eross sections (ie) and (42) have been devised by 


HL 


Tuan, Wright, and Onley [ Ref. re 


D. IDENTIFICATION OF EXCITATION MECHANISMS 


The form factor squared for an inelastic transition is defined by 


do fice) . B? R 
Sal.. or EL ao mott a 


ier the DWBA, the theoretical variation of F4 with Scattering angle can 
be calculated and plotted for each EL and ML at any incident energy. 

In actuality, one needs to consider only L=O through about L=3 since 
higher terms are expected to contribute eer ie amount to the inelastic 
Scattering. Each observed peak of the inelastic eeeere tn can then be 
identified by plotting the experimental form factor variation with angle 
at each incident energy used in the experiment, and comparing with the 


theoretical curves. 
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ITI. EXPERIMENTAL EQUIPMENT AND PROCEDURE 


A. NAVAL POSTGRADUATE SCHOOL LINEAR ACCELERATOR 

The NPS 120 MeV LINAC was described in detail by Barnett and Cunneen 
[ Ref. 13} . It is an electron accelerator similar in design to the Stanford 
University Mark III linear accelerator described by Chodorow et. al. 

[ Ref. 14] . The 30 foot NPS LINAC is powered by three klystrons. Each 
klystron can yield up to 22 megawatts of peak power. The machine is pulsed 
et a rate of 60 pulses per second and each beam pulse is approximately one 
microsecond in duration. The maximum average beam current is approximately 
one microampere. 

Beam deflection magnets and energy defining slits are used to select 
the desired beam energy, which may range from 30 to 105 MeV. The field of 
the beam deflection magnet is measured by a nuclear magnetic resonance 
(NMR) probe. Quadrupole magnets focus the beam to a diameter of approxi- 


mately two millimeters at the target. 


B. SPECTROMETER AND COUNTING SYSTEM 

The NPS LINAC Spectrometer was described in detail by Oberdier [Ref. Tale 
It is a movable, 16-inch ade double~focusing magnetic spectrometer 
capable of measuring electrons scattered at angles from 45° to 150° in 
15 degree stops. 

The field of the spectrometer magnet is varied to measure the energy 
distribution of the scattered electrons. Electrons with energies within 
+1.5 per cent of the spectrometer energy setting ae detected by one of 
ten NE-102 plastic scintillation counters located on the focal plane of 


the spectrometer magnet. Each detector serves as the source of a counting 
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channel. Each counter has a width which allows detection in an energy spread 
of 0.3 per cent of the central spectrometer energy. Two large detectors, 
each of which spans the same field as the ten smaller detectors, serve as 
dual coincidence backing counters. A charged particle passing through a 
front detector and both backing detectors will register as one count ona 
Jorway Model 1836 Dual Channel Scaler located in the control room. The 
counter gating circuit is synchronized with the accelerator gun pulse so 
that those triple coincidences that occur between pulses are not counted. 
The scalers are capable of recording periodic pulses up to 150 MHz. 
However, since the beam is pulsed, the count rate was limited to 50 counts 
per second to avoid significant count rate losses. A Faraday cup is used 
to calibrate the charge collected on the Seconday Emission Monitor (SEM). 
Three pure sla foil targets were used with thicknesses of two mils, 
five mils, and ten mils. The targets were placed in an evacuated chamber 
with 1.0 mil Havar windows positioned to allow the incident electrons to 
enter and the scattered electrons to exit. Since the scattering cross sec- 
tions decrease at the backward angles, five and ten~mil targets were used 
in order to increase the count rates and thereby reduce the effect of sta- 
tistical fluctuations. The two-mil target is recommended Por, fibure and 
vestigations to maintain adequate resolution (the ten-mil 
target gave an elastic peak width of 19.5 per mil while the two-mil target 


gave an elastic peak width of 6.7 per mil.) 


C. DATA ACCUMULATION 
During each experimental "run" the beam energy and scattering angle 
were held constant. The number of counts in each of the ten counting chan- 


nels and the integration capacitance and voltage were recorded automatically 


at each spectrometer energy. The spectrometer energy was decreased in steps of 
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O.1 MeV. The integration capacitance and voltage were chosen to maintain a 
high count rate not to exceed 50 counts per second. The beam current was 
decreased when measuring the elastic peak to avoid an excessive count rate. 
The spectrometer energy was reduced from 100 MeV to about two MeV above the 
elastic peak. Elastic scattering data were taken at 0.1 MeV steps to about 
four MeV below the elastic peak. When the elastic data were accumulated, 
the beam was retuned to give maximum current allowed by count rate restric- 
tions, and the spectrometer energy was set back several MeV to overlap the 
elastic data. In order to insure that the magnet always followed the same 
hysteresis loop, the spectrometer magnet was returned to a field value cor- 
responding to 100 MeV and slowly decreased to the desired energy. It was 
not necessary to consider relative counter efficiencies in this work. The 
KIK code [Ref. 16] automatically corrects the energy of each counter to 


correspond to the true energy being detected. 
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IV. DATA ANALYSIS 


A. DATA REDUCTION 

The data were collected on a Model 33 teletype on paper and punched 
tape simultaneously. The punched tape was transcribed onto computer data 
cards compatible for use on the IBM 360-C. The organization of the data 
cards is described in Appendix A-2. Fortran program DATA LIST was used to 
print out the transcribed data cards for verification against the teletype 


" After verification and correction, if required, the data 


eard Copy. 
cards were processed by fortran program KIK. This computer program assimi- 
lates the ten spectra (one spectrum from each of the ten spectrometer chan- 
nels) into one total spectrum, computes the statistical error, and outputs 
the results in both hard copy and punched card form. The punched cards 
from KIK are used directly by fortran program SPECTRUM to draw the total 
Spectrum as a point plot on the CALCOMP plotters at the NPS Computer Center. 
These spectra are displayed as Figures IV-i through IV-5. 

The spectrum cards from KIK were then processed by fortran program NAW 
[ Ref. 16]. This program best fits an analytic expression to the elastic 
peak data and sone cess under it - the integration is terminated at 
AE = a0, where! is the full-width of the elastic peak measured at half- 
maximum. NAW then makes the appropriate radiation corrections to the area, 
as detailed in Section IV-B. Experience has indicated that the best value 
for parameter L, is L,=3 (see Appendix A-4). 

Program NAW then determines an analytic expression to represent the 
background radiation due to the elastic tail in order to make the requisite 


subtractions. The necessity for these subtractions will be discussed later 


in Section IV-C. Experience of experimenters at Institut Fur Technische 
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a 
en 


Kernphysik at Darmstadt led to the choice of parameter Lo=3 and Lj¢=5 
resulting in the “best"’ representation of the radiation tail. Teme aici 
ation tail and any constant background (measured above the elastic peak) 

are subtracted from the data, and program NAW then determines a best fit 
analytic expression for the inelastic resonance peaks. The technique used 
was to fix the positions of the resonance peaks based on viewing the various 
spectra, select peak widths (full-width at half-maximum) based on other 
published experimental results, and then dynamically adjust these parameters 
until one best fit for all experimental spectra was obtained. Program NAW 
displays these results both graphically and in tabular form. Furthermore, 
it integrates the analytic expression for the inelastic resonances to deter- 
mine the areas under the peaks and displays the ratio of the area of the 
inelastic peak to the area of the elastic peak for each resonance. 

Using the incident electron energies calculated by NAW, fortran program 
MARPHASE uses the Rawitscher-Fischer Partial Wave Code [Ref. 8 | to. -Ccal Cus 
late the elastic scattering parameters; the results are displayed in 
Table IV-1. Since the differential scattering cross sections are propor- 
tional to the area under the resonance peaks one can easily determine the 
inelastic scattering cross sections, or more usefully, the square of the 


2. 


inelastic scattering form factor, F 


2 (do/d5.) inet _ Ai _(de/dr)elas 
Fa (do/d 5) mott ~ Are (do/d.2) mot , 


The inelastic form factors squared may now be graphed as a function of the 
scattering angle and parametized by incident energy. To identify the ex- 
citation involved, one may compare these angular plots with the theoretical 
values as determined by a DWBA code. Fortran program DUKE [Ref. 12] uses 


a DWBA code to calculate Fe for both electric and magnetic transitions. 
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TAREE IV—i 


ELASTIC SCATTERING PARAMETERS FOR wees 


4.30x10-2 


4.91x1072 
3.00x1072 





3.43x1072 


Differential cross sections are in units of pb/sr. 


A transition strength, B, of unity is assumed for all transitions in these 
calculations. While the output labeled as FORM is the theoretical value 
for F* for electric Er aneenen the corresponding output is Fe multiplied 
by the factor V,;/V_, (labeled QVLT on the output) for magnetic transitions. 
Once the theoretical values for F have been obtained the experimental 
values for the transition strength can be calculated by simply averaging 


the values, i.e. 


B= a (Ry / Fa). 


Fortran program BCALC computes B and then multiplies the theoretical values 
of Be by B. These are graphed, and the experimental values of fe are then 
plotted on the “tt graph. This procedure is done with each identified 
resonance for all reasonably possible transitions. The transition type 
identification is made by determining which graph exhibits the best fit of 
the experimental points to the theoretical curves and also yields a reason- 
able value for B. To avoid subjective, and thereby possibly prejudiced, 
evaluations of the best fit, program BCALC computes the normalized quanti- 


ties, VAR, which are compared to determine the best fit - the smaller the 


value of VAR, the better the fit. 


\ 


ZZ 





VAR = => (BiB) 


€ip valves 


These values may be found in Table V-2. 


B. RADIATION CORRECTION CALCULATIONS 

There are processes occurring during elastic scattering that result in 
low energy electrons which constitute the so-called elastic radiation tail. 
Since it is not practicable to integrate the electron spectrum under the 
elastic peak to zero energy, it is conventional to apply three correction 
factors. It is sufficient to correct for such processes as follows: 


(de) . ide) | KKK, . 


iE 


where Ke, Kp, and Ky are, respectively, the Schwinger, Bremsstrahlung, and 
Flectron-Electron Scattering corrections. 
1. Schwinger Correction 

The Schwinger term is by far the most significant of the three 
corrections being discussed. ‘This term corrects for virtual photon emission 
and reabsorption and real photon emission occurring in the field of the tar- 
get nucleus - both of which result in large-angle scattering of the incom- 
ing can i1.¢.5) durinee the scattering event. Schwinger's original 


work [ Refs. Lr 18| gives the relation 


Ke = exe] 42 a 13 | In( 283 sin @/2\- 4 “ree + le (9) sin?e/2) 
AE 12 mc 2 ‘2 2 § 


where ais the Sommerfield fine structure constant, E3 is the most probable 
energy of the outcoming electron, AE is the width of the integration (in 
program NAW AE is set at twice the elastic peak full-width at half-maximum), 


and £(@) is the Spence Function. Many investigations have since improved 
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upon this formulation, notably the work of Tsai [ Ref. 19] and that of 
Meister and Yennie [ Ref. 20] . | 
2. Bremsstrahlung Correction 

This term corrects for the energy radiated via bremsstrahlung during 
small-angle deflections caused by the electron's interactions with the fields 
of nuclei other than the elastic scattering center both before and after the 
elastic scattering event. This energy loss mechanism was first studied by 
Bethe and Heitler [ Ref. 22 and the correction factor for this mechanism 


was given by Hofstadter [Ref. 22| as 


La 


t 
C= exDy = ——— | 
Ks P XoIn 2 Ve : 


where X, is the radiation length and t is the target thickness, both in 


Sf ane. Program NAW uses the formalism of Barber, et. al. [ Ref. 23] to give 


( es Ce a) 1 
K, = : oe eee: |e a oe 
4 ey 2cos ar cCE SED E EE, 3 91n (183exp{-1n z}) 


Pe CeCe eroOn-rleceELon ocattering 


This term, sometimes called the Moller term, is almost negligible. 
The term corrects for electron-electron scattering in the target that occurs 
before or after the elastic scattering event. The effect is to produce low- 
er energy electrons, noes iene Grieinaledirection,  sumilar to the errect 


of the small-angle bremsstrahlung. By Barber's formalism we have 


AZ PRD Maz oo. |) i ] 


C. RADIATION TAIL CALCULATIONS 
When calculating the radiation corrections we were concerned only with 
determining how many elastically scattered electrons were degraded to an 


energy beyond the integration cut-off, AE. The calculation of the radiation 
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tail is Senne with the energy distribution of electrons lost from the 
elastic peak. Program NAW has the option of using either Schiff's method 
[Ref. 24] or the more recent method developed by Ginsberg and Pratt [ Ref. 5] 
The latter was used in this experiment resulting in an approximation to the 
true radiation tail - while the approximation is crude, it does have the 
proper slope. NAW then best fits a fourth-degree polynomial to the cal- 
culated values (output in program NAW as SSW). The resulting expression 
for SSW is then fit to the raw data spectrum after any constant background 
counts have been subtracted. Of the 11 fit options available, the one used 
in this work was 

UNT = A*-SSW + Bt CX, 
where: UNT is the portion of the measured counts/pC attributed to the 
elastic radiation tail; A, B, and C are numerical fitting parameters deter- 
mined by NAW; and X = E-E, (Ep is a constant determined by NAW). The radi- 
ation tail obtained is then subtracted from the inelastic data to yield 
the unfolded spectrum. | 
D. ERROR ANALYSIS 


i Systematic Errors 


In this experiment , the ratio of the inelastic to the elastic cross 
section is computed. Since the inelastic and the elastic data were always 
taken nee the same experimental run, most systematic errors cancelled. 
These would include errors due to variations in solid angle, scattering 
angle, SEM efficiency, target thickness, target purity, and absolute spec- 
trometer efficiency. On several occasions after the inelastic data points 
were collected, the elastic peak points were penaared" tn verify that the 
number of counts per microcoulomb remained constant, within statistical 


variations, during a run. 


ue: 





2. Random Errors 
The statistical fluctuations of the number of electrons scattered 
for a given number of incident electrons obeys Poisson statistics, in which 
the standard deviation is given by S,(E) = y N(E) , where N(E) is the total 
number of counts at scattered energy E. The fractional error is therefore 


given by: 


(N(E) 


Fractional error = 
N(E) 


Program KIK computes the energy, E, for each counter at each spectrometer 
setting, sums the ten spectra to yield one total spectrum, and tabulates 
the total counts per microcoulomb and the associated fractional error. 
This error varied from one percent to three percent. The effect of this 
error on the final form factors is calculated by NAW as follows. Let the 
total counts detected in an energy region encompassing the width of an in- 
elastic peak be N, and the counts attributed to the inelastic peak be Ny 
(thus Ne _ Ny is attibuted to radiation tail and background in this region). 
Then the fractional error in the cross section is given by (Nr Nr. Due 
to the varying strength of each transition and the fact that the total num- 
ber of counts was lower for certain runs, these errors varied greatly in 
Size. They are tabulated in Table V-l. 

Other random errors include changes in the beam current, variation in 
SEM efficiency during a run, energy shifts of the incident beam, leakage 
current in the beam current integrators, error in measurement of the back- 
ground, uncertainties in the position of the peaks, and uncertainties in 
the widths of the peaks. 

Variations in the beam current affect the values obtained for the 
cross sections since count rate losses increase with current. Although 


these losses were initially believed small for count rates up to 50 counts 
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per second, current studies indicate that this is probably not so. The 

magnitude of this error is, as yet, unknown, but may well contribute as 

much as an additional 25 percent error to the experimental values of the 
form factors. 

The possibility of variations of the SEM efficiency during a run was 
not investigated in this experiment. Measurements by Berard and Traverso 
[Ref. 26 | show that the SEM efficiency remains constant within 0.09 percent. 

Variations in line voltage cause fluctuations, not only in beam energy, 
but in spectrometer and beam deflection magnet fields as well. The effect 
is to shift some electrons from one counter to another. However, the elec- 
tronicregulation of the spectrometer and deflection magnets reduce this 
error to a negligible level. Leakage currents in the beam integration cir- 
cuit were measured by Topping [ Ref. 27] and were found to be negligible 
also. 

The background contributes less than 1 percent of the area of the elas- 
tic peak, so that the fractional error due to the background subtraction 
is negligible - particularly when compared to the uncertainty involved in 
the subtraction of the radiation tail. 

The error involved in fitting inelastic peaks to the ane Pe daft icuLe 
to assess. The pee i eacel form factors calculated by program NAW depend 
on the choice of location and width of each of the inelastic peaks. It is 
estimated that these uncertainties contribute a 20 percent error to the 
experimental values of the form factors. 

The wore fractional error for any given calculation of an experimental 
form factor can now be obtained as follows: 


S=/ oS? , 


all sevrces 
et errer 


ay 
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where S; is the fractional error due to the i'th mechanism. Since three 


factors dominate the others, the fractional error can be approximated by 


Ss = Jc.25)2 + (.2)2 + s¢ , where S. is the fractional error due to Poisson 


counting statistics (found in Table V-1). In most cases, the total esti- 


mated error was approximately 40 percent. 
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Figure IV-4 
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V. EXPERIMENTAL RESULTS AND CONCLUSIONS 


A. EXPERIMENTAL RESULTS 

The results of five experimental runs are reported here. Three other 
runs were conducted but are not reported. Two runs conducted at a scatter- 
ing angle of 150° had counts that were so small that any resonances present 
are masked by the statistical fluctuations. One run conducted at an inci-~ 
dent energy of 90 MeV and a scattering angle of 120° using a ten-mil target 
Nadmein Slime lenieresolutdons Ee extrac&eany Structure from the resultant 
spectrum. Prior to presenting the experimental results, a brief description 
of the procedure by which program NAW was used follows. 

The initial positions of the resonances were estimated by viewing the 
experimental spectra, resulting in the values of 7.4, 8.8, 10.6, 13.8 and 
14.2, and 23.0 MeV. Since the fit by NAW for the radiation tail is only 
considered valid for a range of about 20 MeV and the theoretical values for 
the radiation tail are known to become progressively inaccurate with in- 
creasing excitation energy, the 23.0 MeV resonance was not investigated in 
detail (preliminary results indicate that this resonance is probably an 
F2 transition). The remaining four resonances were given to program NAW, 
with widths of 3.0, 2.5, 3.0, and 4.75 MeV respectively, and the program 
was allowed to vary the positions of the last two resonances. Each run 
resulted in somewhat different values for the positions of the resonances, 
and average positions were determined. The averaged values were then en- 
tered into program NAW and, keeping all four resonances fixed in energy, 
the widths were all allowed to be determined by NAW. At this stage, for 
each resonance, the average value of the width for all the runs was calcu- 


lated. Entering these averaged values for the widths into NAW, the 7.4 
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and 8.8 MeV positions were allowed to vary and again the averaging process 
was followed. The new positions were given to NAW, the widths were once 
more allawed to vary, and the results were again averaged. The final re- 
sults of this procedure were resonances at 7.14, 9.12, 10.9, and 13.8/14.2 
MeV with respective widths of 3.13, 2.89, 3.06, and 4.89 MeV. The peak in 
the vicinity of 14 MeV was found to be incident energy dependent. At an 
incident energy of 75 MeV the peak maximum was located at 13.8 MeV, and at 
90 MeV the peak maximum was located at 14.2 MeV. Hereafter this resonance 
will be referred to as the 14 MeV peak. 

Program NAW determined an analytic expression to best fit the data, 
resulting in a smoothed spectrum. In also determines an analytic expres- 
sion for the radiation tail; these functions are displayed graphically in 
Figures V-1 through V-5. The radiation tail and any constant background 
were then subtracted from the smoothed spectrum, and analytic expressions 
for the four inelastic resonances were determined assuming a Breit—-Wigner 
form. The areas under the inelastic resonances were then computed and the 
ratio of each of these areas to the area under the elastic peak was calcu- 
lated along with the fractional error due to counting statistics (Poisson). 
Multiplying A;/A, by the square of the elastic scattering form factor (cal- 
culated by program MARPHASE and displayed in Table IV-1) gives the square 
of the inelastic form factor. The results of these calculations are dis- 
played in Table V-l. 

Using the experimental values of Fe at the different angles and ener- 
gies and the theoretical values as calculated by DUKE as input data, program 
BCALC determines the transition strengths, B, and the values of VAR for 
the various transitions. It also prints out (for subsequent plotting) 
experimental and theoretical values of F* based on the calculated value 


of B. The numerical results are displayed in Table V-2 and the graphs are 


reproduced as Figures V-6 through V-9, 
35 
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B. CONCLUSIONS 

Transition assignments were made as discussed in Section IV-A. The 
7.14 MeV resonance appears to be an M2 transition, but this assignment can- 
not be made conclusively due to lack of scattering data at the background 
angles. Future investigations should concentrate on obtaining data at the 
backscattering angles, since the cross sections for magnetic transitions 
increase at these angles. It may also be useful to obtain data at the small 
forward angles, in order to move away from the peaks in the form factors. 
The heretofore unreported resonance at 9.12 MeV is concluded to be an E2 
transition. The assignment of the 10.9 MeV resonance as an E2 is in agree- 
ment with previous unpublished work at Darmstadt [ Ref. ||. The results of 
our work indicate that the resonance collectively referred to as 14 MeV 
shifts with incident energy. This shift has been explained by Pitthan 
[Ref. 21% who assumed that the broad resonance in the cross section is cal- 
culated from a broad nuclear excitation strength function. This strength 
function varies primarily with excitation energy, but also varies with in- 
cident beam energy and scattering angle. Thus the maximum of the inelastic 
peak may shift when beam energy is changed. The results of this work, 
coupled with the results of the Darmstadt group [4] (conducted at lower 
incident energies) indicate this resonance is an El giant resonance. 

The results described above are recapitulated in Table V-3. All pre- 
viously reported [Ref. 4] resonances appeared in our spectra. The observed 
strengths were not in exact agreement with the previous results, but dis- 
agreement can be explained by observed deficiencies in the counting system. 
Furthermore, two additional resonances were observed: 9.1 MeV, with 
B(E2) = 8,000e2£m4; and 23 MeV which is probably an E2 transition. The 
limitations of the determination of the)elastic radition tail, and insuf—- 


ficient experimental data precluded detailed investigations of the latter 
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resonance. We have demonstrated that giant resonances can be measured at 
the NPS LINAC Facility. The various experimental techniques developed, the 
problem areas identified and rectified where possible, and the computer 
data handling procedures Méioped will enable future investigations of 
giant resonances by the LINAC group to proceed much more rapidly and accu- 
rately. Future investigators should spend a fair oe O time developing 
new techniques for fitting the elastic radition tail, because its effect 


Za 
on the resulting values for Fis paramount. 


TABLE =o 
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Transition 
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APPENDIX A 


COMPUTER PROGRAMS 


Seven FORTRAN programs were used to perform the necessary calculations 
for this work. The computer programs are reproduced on the following pages. 
Each program listing will be prefaced by a brief explanation of what it does, 
and a parameter list with explanatory remarks (program options not directly 
applicable to this and related works will not be listed). Additionally, if 
any special computer control cards are required to execute the program at 


the NPS Computer Facility, a sample FORTRAN deck set-up will be included. 


pom DATA Ls F 

DATA LIST prints the raw data cards in the teletypewriter format used 
by the NPS LINAC Data Acquisition System, segregating by integral data sets, 
so as to aid in verifying the accuracy maintained during the processing from 
punched-tape output to punched-card output. This program requires only one 
parameter card: 

TEST, TEXT —(A4,19A4) 


LESE The first four characters on the first data card. 
Usually a blank followed by three digits. 


TEXT ' Any alphanumeric message to be printed on the first 
page of the data listing. 
LK 
KiK can process up to 499 data sets, each set consisting of the scaler 
readings for each of the ten spectrometer channels of the NPS LINAC Count- 
ing System. KIK compiles the ten separate spectra into one total spectrum 


and calculates the errors, giving the results on both hard copy and punched 
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Carder 


A FORTRAN Deck set-up consists of the requisite control cards, four 


parameter cards, the data deck, and finally two trail cards. 


Ie 


2. 


THETA* 


TEXT Any identifying alphanumeric remarks desired (20A4) 
K, FREO, FREY, DE, TAU; FEA, RH@, THETA (110,7F10.5) 
K=] Program control parameter signifying the start 
of a spectrum. 
FRE@ Upper spectrometer energy limit for the spectrum 
in MeV. 
FREU Lower spectrometer energy limit for the spectrum 
an Mev. 
DF Spectrometer energy step width in MeV. KIK will in- 
terpolate between missing data sets. 
TAU* Counting system dead time/pulse length. 
FEA* Incident accelerator energy (MeV). 
RHP* Target thickness (g/cm2). 
Electron scattering angle (degrees). 


All parameters followed by an asterisk are for identification purposes only 


and are not used in the calculations. They may be left blank if desired. 


She 


aG@lp} ae. a( lo} (1014) 


Spectrometer counting system normalization factors 


relative to counter number 5. The a,'s are defined by: 


a =| ENERGYG) _ | *104 
i | ENERGY (i) 


Input/Output Control 


ipecu. LO, LEOR: INV, FKXADS EKEN CSLe2 PAROS) 

LL=1 

LU, L¢@ The spectrum of counter I specified by LU<I<L@ will be 
punched out onto cards. 

LFOR=0 

INV>0 The order of the input data is reversed; i.e., from 


lower to higher energies. 
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FKXAD>O The input capacitor voltage (full-scale setting) 
will be scaled to units of 107 volts by internal 
multiplication by this factor. For the NPS LINAC 
Facility, FKXAD=100 with the voltage being record- 
ed in ti tlavolts. 


FKEN>O The input energy will be scaled to units of KeV 
by internal multiplication by this factor. For 
the NPS LINAC Facility, FKEN=100 with energies 
recorded in MeV. 
The data cards now follow with two cards per data set. The actual format 
used is quite complex, since the program reads data around the colons and 
reconstructs decimal numbers. The general order of the input data fields 
is as follows: 


Card 1: EG, VFS, E, C, T, V, DUMMY, C@(1)...C@(3) 


Candie? am Gl (4) GO'CLO) > Be A 


E@ Incident accelerator energy 

VES Full-scale voltage setting on the capacitor (mv). 
E spectrometer energy setting (MeV). 

C Capacitance of the collecting cavacitor (uF). 

ae Elapsed time (sec). 

V Final voltage on the capacitor (mv). 


DUMMY = 000000 
CA€1) . Spectrometer counter readings. 
B | Spectrometer backing counter reading. 
A Spectrometer 'accidentals' counter reading. 
De [rade Cardsn 
When a card is encountered with columns 1-7 left blank, that card and 
the immediately following card are interpreted as control cards, with the 
following meanings: 
VFS<O RMIN=VFS*10-3 


=O RMIN=1074 
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Two points differing from each other in energy by 
less than RMIN will be included as one point on 
the spectrum. 


CA(1)>0 Counter I will not be considered in compiling the 
total spectrum. 
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DATA CARDS 


MK tees BKK ESS He 


ERE. 
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PLACE 2 TRAIL CARTS 


my 


He 





ee oF ECLRUM 


SPECTRUM uses the output cards from KIK as its input data, displaying 


the spectrum as a ‘point plot' drawn on the NPS Computer Facility's CALCOMP 


Plotter. 


The graph will be 15 inches long and 9 inches high. SPECTRUM re- 


quires three parameter cards; the first two cards contain alphanumeric iden- 


tifying text in columns 1-48 and the third is as follows: 


NPTS, XSCL, ESCL, EQ@FFST, X@FFST (13273 4810-0) 


NPTS 


XSCL 
ESCL 


EGFF ST 


XOFFST 


Numer of points to be plotted, with the restriction that 
NPTS230m>+ 1; n is®any positive integer. 


Ordinate scale in units/inch (only one significant digit.) 
Abscissa scale in similar fashion. 

Position of the abscissa zero point. Usually for the 
elastic resonance EGFFST=E,, + 2; for the inelastic spectrum 
E@FFST=E,1] which results in the points being displayed at 


excitation energies. 


Position of the ordinate zero point. 
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SAMPLE FORTRAN DECK 
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D. NAW 

NAW uses the output from KIK as its input data. This program calculates 
the ratio of the area under the inelastic resonance peaks to the area under 
the elastic resonance peak. Because of its extreme length (3134 FORTRAN 
cards) the program listing will not be reproduced; documentation may be 
found in reference 15. NAW accepts up to 198 elastic data points, 498 in- 


elastic data points, and 25 inelastic resonances per spectrum. 


1, Text Alphanumeric identifying remarks. (20A4) 
2. EICH, PN,, PN3, PN,, PNc, PN, PNz (7F10.5) 


EICH Spectrometer energy normalization factor to scale all 
input energies to MeV. 


PN Starting values for the various parameters being used to 
fit the elastic spectrum with an analytic expression. They 


are used only when L)>0 on card #7. 


PNo Half-width of a gaussian to fit the high energy side of 


the elastic peak. 


PN2 Half-width of a gaussian to fit the low energy side of the 
elastic peak. 


PN, Matching point for the low energy gaussian and a hyperbolic 
segment. PN, is specified in MeV from the position of the 


zero point and furthermore is a negative number. 


: PN. Position of the zero point. When EN-=0 thegzero poine is 
taken to be the position of the elastic resonance peak. 


PNE>0 Height of the second elastic resonance peak. 
PN7>0 Position of the second elastic resonance peak. 
3. EE, THETA, D, UNT, GINT (5F10.50) 
EE Elastic resonance energy in MeV. 
THETA Electron scattering angle in degrees. 
D Target thickness in allome 


UNT Constant background at the high energy end of the elastic 
resonance spectrum in counts/pC. 
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GINT Lower integration limit for calculating the area under the 
elastic peak. This value is only used when L,=0 on card #7. 
Aesaeerother tames the Lower limit is fixed at Eel minus twice 
the peak width. 


EEL, EEL, GSSW,, GSSW,, DSSW, DWAL (6F10.5) 


EDD eh Cnemeyeliltietor tithing the telasitie peak. 


EEL, Low energy limit for fitting the elastic peak. 
GSSW, Upper limit for calculating the radiation tail. 
GSSW,, Lower limit for calculating the radiation tail. 


DSSW Energy step to be used in calculating the radiation tail. 


DWAL See Ly on card #7. 


Eo? Pas (2H 1025) 
EA Peon energy Limiteton tattane the inelastic spectrum. 
Ea Low energy limit for fitting the inelastic spectrum. 
GR, cri, cr2, cR2, GR3, GR3 (6F10.5) 


Regions where it is expected that only background contributes to 
the spectrum; i.e., between resonances in the inelastic spectrum. 
There must be at least one set GR,, GR, input, with a maximum of 
three such sets. 


Ly, lo, ..-, Log | C27 i) 
Ly Analytic expression used to fit the elastic resonance. ° 
=0 Simple ae sections (AkG3). 


=] A continuous curve consisting of two gaussians and a seg- 
ment of a hyperbola (ARC5). 


eee last Teeresonance wiltebe Eitlas if Ly leand the an- 
elastic resonances will be fit using the Breit-Wigner formula. 
Fach curve will be multiplied by a straight line with equa- 
tion Y=l + ie (12 eal Ome) P; will be specified on card #10; 

eee eeexci tation enerey and let) the position of the 


x s 
maximum. 


pomeeworc idobdest Coolcnel wim eelttedae tt Lat 2 except that 
NAW will determine the values Pj and card #10 will be left 
blank. 


Lo Analytic expression used to represent the background due 
Bomene Pad tation tail 
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=1 BG = P,*SSW 

=2 BG = P1*SSW + Po 

=3 BG = P)*SSW+ Py + P3X 

=4 BG = P, + SSW 

=5 BG = Py + PoX + SSW 

=6 BG = Py + PoX + P4X4 + SSW 

=7 BG = P, + PQA + SSW 

=8 BG = Py + PoA + P3A2 + SSW 

=9 BG = P, + PoB + PB4 + P,X + SSW 

=10 BG = Py + SSW*(P5 + P3X) 

=11 BG = Py + SSW*(P, + P4X + P,X4) 
P. will be determined by NAW. 
SSW = Theoretically calculated radiation tail. 


a eee CLevemecm concur ys wil Bebe intciunal ly edlen— 
lated by NAW. 


= oo =a 1} ® 
A=] (Ed 3) 
B = 1/(E,] - E - DWAL). Should L9=9 be selected, then DWAL must 
be input on card #4. If Lj,#9 then the space for DWAL on card #4 
should be left blank. 
L3 The number of resonances with specified energies. If L.<0 
then the resonances will be kept fixed relative to Eel) 


but E, (1) will be positioned by NAW. 


Ly The number of resonances with energies to be determined by 
NAW. Note that L3 + Ly must not exceed 25. 


Ls The number of line widths to be determined by NAW. 
=Q All line widths will be specified on card #9. 
>0 «6—E.g. =3, means that three line widths are free and in fact 
are Eel ET (238 and E@). ie there areimere than three 
HeS>Olances, Seeger OUCTMG (ijewlil all begiutted wath the 


same line width. 


L¢=0 
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Liy-hy, 


ie He 


Target geometry: 

Transmission —- normal thickness entered on card #3. 
Reflection - normal thickness entered on card #3. 

Effective thickness entered on card #3. 

Output control: 

Maximum output. 

Minimal output - no tables and no intermediate calculations. 
Limited output - no tables. 

PLOceim.e Control: 

Results are not plotted. 

Results are plotted. 

Repetitions of the NAW procedure. 

No repetitions are performed. 

In this case it is possible to evaluate the data with new 
parameters. The new parameters will be entered on additional 
cards placed after the data. 


Inelastic spectrum fits iteration control. Only override 
to the default options need be entered. 


Maximum number of iterations; default is Ly9=30. 


Minimum number of iterations; default is 11375: 


The radiation tail will be calculated using Ginsberg's method. 
The radiation tail will be calculated using Schiff's method. 


Radiation tail fit control. Only override options need be 
entered. 


Analytic representation of the theoretically calculated 


Yadiation tail: default as Ligul- 


Ze 6 4 
Po of PoA a PAA BA oF PoA 


g) 2 
Py + PoA +} PoA + PLE = Pok 


SSW 


SSW 
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cae 


Bb 


SC. 


fhe 


10. 


~ a 272 

=3 SSW = Py + PoA + Pak + PARES + PoE3 

an Z =) 2 aa 

=4 SSW = PLA + PoE + PE + Py A + PrE 
5 SSW = P, + PoU + PU + PAU + PeU 


<O 


Li7 


Lig 
lig ale 


Lo9=0 


Ly 1>0 


ae will be determined internally by NAW. 
A = TRGB) al E) 


mies (b 


= GSSW, + (GSSW, - GSSW,)/Ly. 


U 


a 


The radiation tail parameters P, and E, will be specified 
Onmeard 15), 


The number of parameters for the radiation tail fit 
(3 $1,555); default is L)5=4. 


see Li6; default is Li g=2- 


Error ellipsoid control values will not be calculated near 
Eheaf tt. 


The calculation of the radiation tail will be modified by 
the factors entered on card #16. 


L >0 Inelastic spectrum results will be punched out onto cards. 


Ze 


E.(1), 


@), soy MCs 102 (7F10.5,110) 


Resonance excitation energies in MeV; enter those with 
variable position first. For more than seven resonances 
set LE>0 and then continue listing the additional reso- 
nances on successive cards. 


E,(8) .-- He(14), LE 


2G) peo EC) 5 Hs 
LCD) gual ee 
Wal) NM) wc ee V7) oe LDR (7F10.5,110) 


P(1) 


Resonance line widths. For more than seven resonances set 
LBR>O and continue on successive cards as on card #8. If 
some or all V(I) are zero, then the line width used will be 
the dispersion of the spectrometer counting system. 


Ore 2 Cs Le CAMO MILI G)), 


The P; factors used when L)=2 on card #7. Set LP>O for 
more than seven factors and continue as for card #8. 
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ng hs 


2 


ee 


tere 


PraAg KAD, hilo, Ky, G, CG, I (OnPOns)) 

ib Atomic number. 

A Mass number. 

XRAD Radiation length in g/cm 

RMS Root mean square nuclear radius in fm. Used only when N=2 
on card #14, otherwise left blank. 

R, G When N=0O or N=1 on card #14, R and G are interpreted as 
the parameters in the Helm expression for the form factor. 
When N=3, R=sCRT and G=TRT. When N=2 or N=4, R and G are 
left blank. 

C, T The fermi charge distribution parameters. Since these val- 
ues are used as the parameters for the trapezoidal charge 
distribution potential used in the radiation tail calcula- 
tions, the values must be entered even when N#4. 

XMAS, H (2F10.5) 

XMAS Nuclear mass equivalent in MeV. 

H Isotope percentage (100% = 1.00). 

3091s ame Oommen ce GE10e=) 
The A... E are the coefficients for the calculation of the 
elastic form factor using a power series expansion onq : 
a See 4 6 8 
fi ees ge oe qe tla sag 
When A=0 the elastic form factor will be calculated using 
the Born Approximation. 
A, N CPS 24x13) 

A - Mass number. 

N Selection of the nuclear model to be used in the Born 
Approximation to calculate the form factor: 

=0,1 F= 3(4 1 (qR) /qR) exp (-G4q2/2) Helm. 

=2 F = 1 -q” <R*>/6 + q‘<R°>“/120 + q°<R2>°/5040 Rethe. 

=3 Fe= (1 - CRT* q~) exp (- TRT* q7) Shell model. 

=4 The form factor will be calculated using the Fermi model. 
EL eee. 4. Meee. Ey (6810-2) 


ibe Cakeders reatisenty aie). <0: 
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Gi: 


SSW 


PI-P5 Radiation tail parameters, see Lig for the expressions. 
Ey This value is read only when Lj¢=-5. 

XSSW, XM, XB, LSSW GEIO.55 23) 
This card is read only if Lo)>0. 
XSSW Factor multiplying the area under the elastic peak. 
XM Factor multiplying the Moller term. 
XB Factor multiplying the Bremsstrahlung term. 


The equation for the radiation tail is now given by: 
(Max IL) 
2 
F On 
LSSW20 The calculated radiation tail will be punched onto cards. 


= ( + XM(Moller term) + XB(Brehm. term) )* (XSSW)A, 


The data cards now follow. The elastic and inelastic data must 
be separated by a card having the data energy step width entered 
beginning in column 72. Repeated evaluations of the same set of 
data. (ite, L179) must have the following set of cards placed 

after the data cards - one set for each repetition to be performed: 


1) Ib beer see card #7 
2) Vo My see card #5 
BE. Ce. BN) see card #8 
4) V(l) ... VN) see card #9 
SIO EGL) FR ea), see card #10 
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E. MARPHASE 
MARPHASE uses the Rawitscher-Fischer phase shift code YALERF developed 
at YALE [ Ref. 8] to calculate the elastic electron scattering cross sections, 


Mott cross sections, and form factors. A parameter set consists of two 


cards. 
lee A een 2 PS... MODEL C5F 10.5110) 
Z Atomic number. 
A Mass number. 
P3=0 
MODEL=l1 Fermi model of the nucleus is assumed with 
a 
[+ exp[r~¢/(t/4.4)} 
In this case Pl=c and P2=t, both in fermis. 
=2 The Shell model of the nucleus is assumed with 
oye ale 
o(r)= [i+ a¢Q)'| expf-(4)) . 
In this case Pl=a@a and P2=A. 
ee ieee eT. RO, LAG (SF10.5) 
E@ Electron beam incident energy in MeV. 
Tl Starting scattering angle in degrees. 
dW Final scattering angle in degrees. 
TH Angular increment in degrees. 
RO The fitting radius pertaining to matching the Nuclear 


and the Coulomb wave functions. 
LAG The number of phase shifts. 
3. Another card #2 placed next will cause the program to use the 
same information for card #1 to compute the new combination of incident 


energy and scattering angles. A blank card is interpreted to mean that 


a3 





a computation is now desired for a different set of card #1 parameters and 


is followed by a new parameter set, usw. 
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Bev UKE 


DUKE caleulates the theoretical cross sections for both) electric ‘and 


magnetic giant resonance transitions. Either the DWBA or the PWBA may be 


selected. A parameter set consists of four cards. A subsequent parameter 


set consists of only the last three cards. A blank card signifies the end 


of the desired calculations. 


Z, AX, Zl, El, E2, ETRAN, BORN, MORE 7 Giese alee) 
Z Atomic number. 

AX Parameter c for the Fermi charge distribution. 

ZA Parameter t for the Fermi charge distribution. 

El starting energy in MeV. 


EZ Final energy in MeV. 

ETRAN Excitation energy in MeV. Set ETRAN=0 for a pure coulomb 
transition and set ETRAN=F1-E2 if the transverse component 
of the cross section should be considered. 

BORN=0 Distorted Wave Born Approximation. 

=] Plane Wave Born Approximation. 


M@RE=O0 Electric transition. 


=| Magnetic transition. 


. L, KMAX, RADVAR, SKIVAR, ASG, ZKMAX, ETRINT, LM@D (2110,5F10.5,110) 


i Multipolarity of the transition. 
KMAX>O Number of partial waves. 

=Q KMAX is calculated internally: KMAX =9+ (E1-30)/10. 
RADVAR Position of the transition distribution normalized to c. 
OKIVAR Width of the transition distribution normalized to t. 
ASG=0 Complete output is printed. 


=leeo Limited output Wo printed. 


90 





ZKMAX This modifies KMAX, as computed internally, to be 
KMAX = KMAX + ZKMAX. 


ETRINT Parameter used to determine the integration steps. ETRINT = 
35.0 by default; increasing ETRINT decreases the step width, 
increases the accuracy and increases the execution time. 


LM@D=0 Transition distribution function will be approximated by 
the derivative of the ground state fermi distribution func- 
tion. 


=] Transition distribution function will be approximated by the 
ground state fermi distribution function. 


Hn ee Oo li NSE NST 3 tot | (2F10.5,4110) 
This card is best given as a blank card, DUKE then determines the 


integration parameters using the internal SUBPROGRAM INDEX; see 
below for assumed values. 


HH1 Integration stepwidth inside the potential in fermis. 
HH2 Integration stepwidth outside the potential in dimensionless 
units. 


NSTl Length of the internal region in units of AX. 


NST2 Length of the external region in dimensionless units over 
ETRINT. 


NST3 The values of the radial wave functions will be computed in 
three=recions delimited by. NSTI aul + NST2, HH2 + i°NST3.5 oi 
wheelies =i e2ee 2.ceeelo . le 


When this card is left blank then: 


oleae Ow fits . 

HH2 = 0.3he/ J E1+E2’ 

NST1 = 3AX/HH1 

NST2 = (30. + (ETRAN/E1):+ETRINT) /HH2 


Seed Mean of the integrated radial wave function within one wave- 
TST) length of the? oscillation with the Wargest wavelength. IST1 = 6, 
and NST3 = wavelength/ (HH2-IST1) 


THEMN, DELTH, N@, STANZ, NZDG PET CE gl lly WO pikes a ya) 
The cross section for a parameter set will be calculated for all 
angles 9 = THEMN + i:DELTH, i=2, 1, ...,N@. KMF, STANZ, and KZDG 


are test parameters and have no effect when left blank. If this 
eiramismlefteblank then THEMN = 30, DELTH > 15, and NO = 10. 
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G. BCALC 

BCALC is used to calculate the transition strength, B. While written 
specifically for use with 75 MeV and 90 MeV incident energies, it is easily 
changed to accommodate other values. BCALC uses the following data cards: 
one card with the number of 75 MeV experimental points; a card for each such 
point containing the scattering angle in degrees and the experimental value 
eit Ba. FORMAT (13,7X,E10.4); a similar set-up for the 90 MeV experimental 
points; and lastly, the theoretical values for scattering angle and F?, 
For each transition type to be evaluated a set of theoretical values must 
be supplied - 75 MeV values followed by 90 MeV values. The output from 
BCALC is a listing of all experimental values with the corresponding the- 
oretical values printed next to them, the computed value for B, the com- 


puted value of the quantity VAR (see Section IV-A), and finally a listing 


of the theoretical value of F4 multiplied by B. 
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